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Abstract: We generate charged black brane solutions in D— dimensions in a theory of 
gravity coupled to a dilaton and an antisymmetric form, by using a Harrison-type trans- 
formation. The seed vacuum solutions that we use correspond to uplifted Kaluza-Klein 
black strings and black holes in [D — p) -dimensions. A generalization of the Marolf-Mann 
quasilocal formalism to the Kaluza-Klein theory is also presented, the global charges of 
the black objects being computed in this way. We argue that the thermodynamics of the 
charged solutions can be derived from that of the vacuum configurations. Our results show 
that all charged Kaluza-Klein solutions constructed by means of Harrison transformations 
are thermodynamically unstable in a grand canonical ensemble. The general formalism is 
applied to the case of nonuniform black strings and caged black hole solutions in = 5, 6 
Einstein-Maxwell-dilaton gravity, whose geometrical properties and thermodynamics are 
discussed. We argue that the topology changing transition scenario, which was previously 
proposed in the vacuum case, also holds in this case. Spinning generalizations of the 
charged black strings are constructed in six dimensions in the slowly rotating limit. We 
find that the gyromagnetic ratio of these solutions possesses a nontrivial dependence on 
the nonuniformity parameter. 
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1. Introduction 

Finding new solutions of the theory of general relativity by using solution generating tech- 
niques is a subject of long standing interest. One of the main approaches consists in the use 
of hidden symmetries (dualities) arising in dimensionally reduced theories. This technique 
was initially used in four-dimensional general relativity coupled to a Maxwell field Q and 
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was generalized in recent years to higher dimensional gravity coupled to various matter 
fields, a great variety of important solutions being generated in this way. 

This approach is particularly relevant in the case of Kaluza-Klein (KK) solutions with 
a nontrivial dependence on the extra-dimension. The configurations of interest here are the 
static nonuniform black strings (NUBS) and the KK localized black holes (BHs) (thus we 
shall not consider KK solutions with bubbles Q, |p, Q or multicenter black objects ||5|). 
The numerical results in the literature reveal in this case the existence of a complicated 
phase structure, the caged BH and the NUBS branches merging at a topology changing 
transition @. Most of the results in the literature on KK solutions with a dependence 
of the extra-dimension concern the case of vacuum configurations. Even in this case, 
there are no explicit solutions that one could write down in closed analytical form (both 
BHs and NUBSs cannot be found by using a Weyl ansatz [^] since the number of abelian 
Killing vectors is too small for them to belong to a generalised Weyl class). However, one 
can analyse their properties by using a combination of analytical and numerical methods, 
which is enough for most purposes. 

The more realistic case of localized BHs and NUBS solutions with matter fields is rel- 
atively scarcely explored. Since it is extremely difficult to solve the field equations for such 
configurations, it would be very helpful if one could generate new solutions from known 
ones on the basis of internal symmetries of the field equations. The main purpose of this 
work is to present such a construction for charged black brane configurations in a theory 
of gravity coupled to a dilaton and a {p + l)-form and formulated in D-dimensions. These 
solutions are generated by using a Harrison-type transformation originally proposed in Q, 
applied to the known KK vacuum nonuniform black strings and black holes m D—p dimen- 
sions. Different from other cases previously discussed in the literature |p, this approach 
produces solutions for an arbitrary coupling of the dilaton field to the antisymmetric form. 

To analyze the thermodynamic properties of such non-uniform brane solutions one has 
to be able to compute their total action and define the various conserved quantities. For 
asymptotically flat spacetimes there are several ways to compute the conserved charges 
and action. One particularly interesting approach is the so-called counterterm method, 
which was inspired by the AdS/CFT conjecture in string theory. In this approach one sup- 
plements the action by including suitable boundary counterterms, which are functionals 
of curvature invariants of the induced metric on the boundary. By choosing appropriate 
counterterms that cancel the divergencies in the total action, one further deflnes a bound- 
ary stress-tensor by varying the on-shell action with respect to the boundary metric. Using 
the boundary stress-tensor one can then construct a conserved charge for any asymptotic 
Killing vector. However, unlike the AdS case, here one does not have a universal prescrip- 
tion for deflning the boundary counterterms and their form seems to depend heavily on the 
boundary topology. A novel counterterm for asymptotically flat spaces has been proposed 
recently by Mann and Marolf Q. In contrast to the previous counterterms, by its definition 
the Mann-Marolf (MM) counterterm assigns zero action to the Minkowski spacetime in any 
dimension, while the conserved quantities are constructed essentially using the electric part 
of the Weyl conformal tensor. In this work, we also propose a generalization of the MM 
counterterm to extended objects, like the p-brane spacetimes. These spaces are not asymp- 
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totically flat but ^transverse asymptotically flat\ Leaving a more general analysis for future 
work, here we show that the total action and conserved charges for p-branes can be defined 
in terms of the electric part of the Weyl tensor and the final results are consistent with 
those previously derived (by other means) in the literature pO| , 11, 12, 13]. As a check we 



also compute the action and conserved charges using different choices of counterterms for 
asymptotically flat spaces and we obtain similar results. One should note here that, unlike 
the case of the MM counterterm, these other counterterm choices are severely limited as in 
higher than eight dimensions and in certain coordinates they cannot regularize the action 



even for flat space |14f| . 



Once one computes the conserved charges and total action, one is ready to perform 
the standard thermodynamic analysis of the generated solutions. In particular, we proved 
that the thermodynamic properties of the charged solutions generated using the Harrison 
transformation can be simply derived from those of the vacuum seed configurations. The 
neutral and charged variables are related in a simple way. Moreover, we find that the 
thermodynamics of the system is affected by the value of the dilaton coupling constant. 
Our results show that all D = 5,6 charged KK solutions constructed in this way are 
thermodynamically unstable in a grand canonical ensemble. However, when considering 
instead a canonical ensemble, one finds stable solutions for large enough values of the 
electric charge. 

The outline of the paper is as follows. We begin with a presentation of the general 
ansatz and we give the necessary details regarding the Harrison transformation in the next 
section. In Section 3, we consider the computation of the global charges of black branes, 
proposing a generalization of the MM quasilocal formalism to the case of the KK theory. In 
Section 4 we review the basic properties of the caged BHs and NUBS vacuum configurations 
of the KK theory that are used as seed solutions in the Harrison transformation. Section 
5 contains a general discussion of the charged black branes in KK theory. A more detailed 
analysis is presented there for D = 5, 6 NUBS and caged BHs in Einstein-Maxwell-dilaton 
theory. The static, charged NUBSs generated in this way are used in Section 6 to construct 
D = 6 spinning configurations. The results there are valid in the limit of slow rotation. 
We give our conclusions and remarks in the final section. The Appendix A contains a 
brief discussion of some technical aspects involved in the numerical construction of the 
caged BHs and NUBS vacuum solutions, including an independent confirmation for the 
existence of D = 6 caged BHs. In Appendix B we present an exact solution describing 
slowly rotating, charged uniform black strings in EM theory in D dimensions, and compute 
the gyromagnetic ratio of these configurations. 

2. General formalism 

2.1 The action and field equations 

We consider the following action 

I = ^ [ d^'x^ (r- -iV<pf - -^^F^,^ f d^'-^xV^K, 

IQttGJm ^ \ 2^ ' 2(p + 2)! p+V SttG Jqm 

(2.1) 
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with (/) a dilaton field and a (p+2)-form field strength with Fp+2 = dAp+i, where Ap+i 
is the corresponding (p + l)-form gauge field. The free parameter a governs the strength 
of the coupling of the dilaton to the form field, hij is the induced metric on the boundary 
dM and Kij is the extrinsic curvature of this boundary, with K = Kijh^K For suitable 
values of -D,p, a, the action ( |2.1D corresponds to truncations of various supergravity models 
and also of the bosonic part of the low energy action of type IIA and IIB string theory. 
The corresponding field equations are 

d.iV^gg^^dA) = -J^_^~"^'Pp^2, (2.2) 
d, (e-2«^V=^F;;--^+^) = 0, 

where the total stress-energy tensor has contributions from the dilaton and the (p + 2)-field 
strength: 

Although we shall consider in this paper only solutions with an 'electric' potential ^p+i, 
magnetically charged solutions can be obtained by using the generalized 'electric-magnetic 
duality' by taking: 

a(l)^-a(l), p^D-p-4, Fp+2 ^ Fd-p-2 = e^""^ * Fp+2, (2.3) 
where * stands for the Hodge dual, the geometry being not affected by this transformation. 
2.2 The Harrison transformation 

Following the method outlined in Q (see also ) , we suppose that the space-time has (at 
least) p+1 commuting Killing vectors, one of them being timelike. Then the D— dimensional 
line element can be written in the following form: 

ds^ = glf{x)dy'dy^ + g^^l{x)dx''dx^ , (2.4) 

where i, j = 1, . . . ,p+l, while a, f3 = 1, . . . , D—{p+l); also note that go = det{g^^^ (x)) < 0. 
The above line element is assumed to be a solution of the vacuum Einstein field equations. 

Upon performing the KK reduction with respect to the y*— directions, the dimension- 
ally reduced action corresponds to a non-linear (D — p — l)-dimensional u-model with the 
target space SL{p + l,R)/SO{p + 1) x SL{2, R)/ S0{2) x R. These symmetries imply 
the existence of a Harrison transformation nontrivially acting on the spacetime variables 
and matter fields. Therefore one may generate nontrivial solutions of the equations ( |2.2| ) 
starting with known vacuum configurations. A detailed description of this procedure and 
the explicit form of the Harrison transformation at the cr- model level is given in ref . 0] . 

As a result, one finds a nonvacuum solution of the equations (|2.2| ), with the line-element 

ds^ = n~-^<^-'P-^'^gf^{x)dy'dy^ + 02"(P+i)^7j°J(x)(ix°(ix'^, (2.5) 
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the dilaton field and the only nonvanishing component of the {p + l)-form being 



= 00 - 2aa(Z) - 2) log 17, (2.6) 

Aty,..y, = V'^{D-2)ae'"^°\go\n-^U, (2.7) 

where 

a = [2a\D-2) + {p+l){D-p-3)]-\ and n = {I - U^)-\l - U^\go\), (2.8) 

U being the parameter of the Harrison transformation, with < U < 1. Here (po is a real 
constant corresponding to the asymptotic value of the scalar field; in what follows we shall 
set 00 = 0, without any loss of generality. 

In this way, one finds a family of solutions of the Einstein-dilaton-antisymmetric form 
equations, parameterized in terms of U. As we shall see later, the parameter U physically 
corresponds to the electrostatic potential difference between the horizon and infinity. One 
should also remark at this point that the new line element (^]^) preserves the symmetries 
of the seed metric (|2~ 



2.3 The metric ansatz and asymptotics 

The procedure described above works for any vacuum solution of the Einstein equation 



that can be written in the form (2^), regardless of its asymptotic expression. However, a 
(p+1)— form gauge potential naturally couples to the world volume of a p— brane. We recall 
that the black p— brane solutions are rich dynamical systems, exhibiting new behaviours 
that have no analogues in the black hole case |16|. 

In this work we consider solutions approaching at infinity the Minkowski spacetime 
M^^P^^ times a Ricci flat space x S^. The metric ansatz that we consider is 

ds^ = gtt{r, z)dt^ + grr{r, z)dr^ + g^zir, z)dz^ + g^^{r, z)du;l)_p_^ + gyy{r, z) (2- 

with r the radial coordinate in M^~p~^ spacetime and duj"^ is the unit metric on S"^. Each 
y -direction is compactified on a circle of length Li (we denote Vp = nf=i ^i)- ^ ^ 
periodic coordinate with period L that parametrizes S^, while t is the time coordinate. For 
any topology of the event horizon, the metric functions possess the following asymptotic 
behavior, with a power law decay^ 

gtt^ -1 + ^U^, gyy^l + ^J^, 9zz^l + ^T^, ~ r2(l + O(l/r)),(2.10) 

Q, Cy and Cz being real constants, which as discussed in the next section, encode the global 
charges of the solutions. 



^The situation is more complicated in D—p = 5 dimensions for black objects with an explicit dependence 
on the 2— direction. In these cases there are logarithmic terms in the asymptotic expansion of the metric 
functions. For g^uj, this occurs already at the next to leading order term, g„„ ~ r^(l + ulogr/r), with u a 
term that depends on Ct and c^. However, the first two terms in the expression of gtt,gzz and gyy are still 



given by (2.1C). More generally one can have different constants Cy for each Y-direction; however this will 



not be the case for the solutions considered in this paper. 
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The generic asymptotic expression of the electric {p + l)-form is given by 

Qe 
.D-p-4 



■^ty,..y,=^-J^ + ..., (2.11) 



while the dilaton field decays as 

^ = -;i^.+---' (2-12) 

with qg, real constants. 

3. Conserved charges and the quasilocal formahsm 
3.1 The Mann-Marolf counterterm 

In general, the conserved quantities of black P-brane spacetimes are encoded in the asymp- 



totic expansions of the metric functions, as given for instance in (2.10). Their computation 
is a problem discussed by many authors, starting with [|lC| ] (see also |11, |l^, |l^), mainly 
by making use of a reference background, whose obvious choice in the case of a general 
P-brane is A4^^^ x T^, where is the {D — i-')-dimensional Minkowski spacetime 

(to simplify the relations, we denote in this Section P = p + 1). 

An alternative procedure motivated by the AdS/CFT conjecture and referred to as 
the 'counterterm method', which avoids the difficulties encountered in the background 
subtraction method, makes use of the quasilocal stress-tensor of Brown and York [^]. In 
this approach one supplements the action by including a suitable boundary counterterm: 

Ict = ^ [ d^'-^xV^k. (3.1) 



dM 



As in the better known AdS case, the counterterms K are functionals only of curvature 
invariants of the induced metric on the boundary and so they do not affect the equations 
of motion. Conserved quantities can be constructed using a boundary stress-tensor Tij 
obtained via the algorithm of Brown and York, that is by taking the variation of the total 
action (including the counterterms), with respect to the boundary metric hij. Provided 
that the boundary geometry has an isometry generated by a Killing vector a conserved 
charge 

Q^ = £ d^'-^S' eTij (3.2) 

can be associated with a closed surface S. Physically, this means that a collection of 
observers on the boundary with the induced metric hij measure the same value of Q^. 
Although the counterterm method gives results that are equivalent to those obtained using 



the background subtraction method [14, [l^], we employ it here because it appears to 
be a more general technique than background subtraction, and moreover, it is interesting 
to explore the range of problems to which it applies, in particular for configurations with 
a dependence on the compact extra-dimensions. 



^When one can apply the background subtraction method. 
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Recently, a particularly interesting counterterm for asymptotically locally flat space- 
times has been put forward by Mann and Marolf g (see also [|l8|, H, |2|). This 



counterterm is taken to be the trace K oi a, symmetric tensor Kij , which in general dimen- 
sions^ is defined implicitly in terms of the Ricci tensor TZij of the induced metric on the 
boundary via the relation: 

= k,,k -k'^kkj. (3.3) 

In solving this equation one picks a solution k^j that asymptotes to the extrinsic curvature 
of the boundary of the Minkowski space as the manifold boundary dM is pushed to infinity. 
One should note here that the boundary data (9A^, h) is not uniquely defined by the bulk 
spacetime {M.,g) but it also depends on the limiting procedure used to define the Gibbons- 
Hawking term. As shown in ||2^, if one employes the so-called 'cylindrical cut-off', then 
the renormalized action leads to a boundary stress-tensor whose leading-order expression 
in D > 4 dimensions involves the electric part of the Weyl tensor: 

where Eij is the pull-back to the boundary dM of the contraction of the bulk Weyl tensor 
Cfiupr with the induced metric h^i,. More precisely, introducing the unit normal vector 
to the boundary dAi then the electric part of the bulk Weyl tensor is defined by: 

E^u = C fipvTfi^n = ~C ^pi,rh^ , (3.5) 

while Eij is simply the pull-back to the boundary of the above tensor. In the 'cylindrical 
cut-off', which is more suitable for practical computations, the asymptotic form of the 
metric at spatial infinity takes the form: 

ds^ = - (1 + 0(r-(^-3))) dt" + (1 + 0(r"(^-3))) dr^ + (a;,j + ©(r'^^-^))) dO' dO^ 
+0{r-^^-^^)de^dt, (3.6) 

where lojj and 9^ are the metric and coordinates on the unit (D — 2)-sphere. 

Our aim here is to generalize these results and compute the action and conserved 
charges for the general P-branes considered in the previous Section. These spacetimes are 
not asymptotically flat and they belong to the class of the so-called ^transverse asymptoti- 
cally flaf spacetimes. In general, one can compute the conserved mass, angular momenta 
and brane tensions for these spaces by using background subtraction procedures if one 
takes the reference background of the form A4^~^ x [O, 12, |l3|. Since the background 



spacetime is flat, one expects to be able to extend the construction of the Mann-Marolf 



counterterm to apply for this case as well and, following the steps in [^, we shall proceed 
to do so. 



^In four dimensions its explicit form has been recently derived by Visser pl| ]. 
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Generalizing the ansatz in (3/7) to apply for P-branes in D-dimensions, one takes the 
asymptotic form of the metric at spatial infinity of the following form: 



a=l 

+ 0{r-^''-''-^^)dy-dy' + (u^u + 0(r-(^-^-3))) dO' dO' + O(r-(^-^-4))d0^dt, 

(3.7) 

where w/j and 0^ are the metric and coordinates on the unit (D — P — 2)-sphere. Here r 
plays the role of a 'radial' coordinate such that spacelike infinity is reached in the r — > oo 
limit and the symbols 0{r^^^^^^^^) refer to terms that fall-off at least as fast as 7--(^~^~3) 
in the r ^ oo limit. 

Note that the extrinsic curvature Kij satisfies the Gauss-Codazzi relation: 

T^ij = Rikjih^^ + KijK — K^Kjm, (3.8) 

where Rikji is the pull-back to the boundary dM of the bulk Riemann tensor. The leading 
terms of Kij and Kij are given by the extrinsic curvature of the standard cylinder of radius 
r in the reference background M.^~^ x T^: 

Kij ~ kij = roj,, + 0(r-(^-^-^)), (3.9) 

where ujij is the pull-back to dAi of the round metric loij on the unit sphere S^~^~'^. 

As shown in |l^, in order to compute the difference Ajj = Kij — Kij to first order 
one has to solve the linearized Gauss-Codazzi relation: 

Rikjih^^ = -Akih'^^kij - Aijk + Auk^j + Ajikl . (3.10) 

Recall that hij is the induced metric on a hypersurface r = const, with normal n* such that 
9ij = hij + riiTij. It is convenient at this stage to further define ^ij = Kj + UiUj — Y^^^^ yfUj, 
where is a future-directed timelike vector on dM, associated with a foliation = {t = 
const.} of the spacetime Ai, respectively (y")* is an outward pointing spacelike vector 
associated with a foliation S^a = {y-^ = const.} of M. Then kij = ^ + ©(r'^^-^-^)), 
while k = + 0{r~^^~^~^^) and replacing these relations in ( 3.10D one obtains: 

+Z. ^ 



a=l 

p 

+^{hij + UiUj - yty^) = -Rikjih''', (3.11) 

a=l 

where we denote Aj = AijU^ , A = AijU^u^ , respectively Af = Aijy°-^ and A"'' = Aijy°-^y^K 
Quite generally, we shall consider matter sources for which the stress-energy tensor at 
infinity falls-off fast enough such that the Riemann and Weyl tensors agree to leading 
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order at spacelike infinity (this is the case of the solutions in this work). Therefore, one 
can replace in the right-hand side of the above equation the Riemann tensor by the Weyl 
tensor such that the right-hand side becomes equal to the electric part of the Weyl tensor 



Eij. 



Contracting now ( |3.11| ) with one obtains: 



— ; — - V + ^ r = ^^^-^ ^^-^^^ 

a=l 



and further contractions with u^, respectively y^^ leads to: 

D-P-2 ^ D-P-2 



^ = r. r. ^ Ei,v}v?, A" = ^ ^ Ei.uY' (3.13) 



respectively. Similarly, contracting (3.11) successively with y"* and y^ leads to: 



Finally, taking the trace of ( |3.11| ) with /i*-' one obtains: 

A = ^ ^"=1^ . (3.15) 
D-P-3 ^ ' 

One is now ready to compute the boundary action (for vacuum spacetimes this is the 
total action): 



P 



= -^L """^^ ^D - P - 2)(0 - P - 3) ^'^'"'"^ - g 

Now, in order to compute the conserved charges we evaluate the boundary stress-tensor 
using the formula (see eq. (4.2) of (P): 

T^J = -^i\J-h,,A), (3.17) 

from which the conserved charges associated with a boundary Killing vector ^ are found by 
direct integration as in ( |3.2D . The conserved mass is generally defined using an asymptotic 
time translation Killing vector ^ = d/dt, while the angular momenta are defined using 
angular Killing vectors ^ = d/dip. If Y^t is a spatial slice and is its intersection 

with the boundary dM at spatial infinity, then one integrates (3^) on over the 

brane directions and the 5^-^-2 gpj-^grg infinity to compute the mass and angular 
momenta. Therefore, to compute the conserved mass and angular momenta one basically 
has to evaluate the quantity: 

STiGTijU^ = -{Ai - Aui) (3.18) 



(D-p-2)(zj-p-3) - ^ - 2)^^.-^' + E Eimv^'iy^'^^^ - ^""yf)] ■ 
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where we made use of the formulae ( 3.12| ) and ( |3.15| ). As an example, the conserved mass 
will be given by: 

M = ^ / d''yd°-''-H^T,fiu> (3.19) 

where ^ = d/dt and Vp = f d^y is the volume along the brane directions. Similarly, the 
angular momentum will be computed using: 

J=^j d''yd''~''-^^Ti,eui (3.20) 

where now ^ = d/d(p. 

The gravitational tension tensor can be defined similarly using the asymptotic spatial 
translation Killing vectors d/dy"" and the integral is now taken over all the spatial directions, 
excluding however the brane direction [23|. To compute the tension tensor density one has 
to evaluate the quantity: 



87rGT,,y-^ = -(Af - Ay,") 



r 



■[{D-P- 2)E,,y-^ + Eirny'^'in^n^ - y'^^yf) 



{D-P-2){D-P-3) 
p 

-EUu'u"'-Y^y^ly'>--)y-]. (3.21) 
b=l 

In consequence, the tension tensor density (defined as in ||2^) is given by: 

T'^' = I d''-''-H^T,,y-'y'^ (3.22) 

' d^'^^'eVJI ,^_,_l^_^_,, [{D-P- ?.)E,,y-y'^ 



SttG J ^'^{D-P-2){D-P-3) 

p 

c=l 

Note that T"'^ is a symmetric tensor. 

Suppose now that, regardless of the topology of the event horizon (if present), some of 
the metric functions appearing in (3.7) have the following asymptotic form^: 

gtt ~ -1 + ^£,_p_3, fl't/^y'" ~ 1 + ^jj_p_^, gyayb ~ ^£,_p_3, (3.23) 



Note that the asymptotics (3.23) are more general than (2.10), since extra-diagonal components in 
the j/'^-directions are allowed. This would cover also the case of KK solutions with an arbitrary torus 
compactifications, so far considered in perturbation theory only [M. 
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where Cq and dab are real constants, with a = 1..P. 

Then it can be checked that to leading order we have the following relations: 

_ {D-P-2){D-P-3) ct j^ya_ {D-P-2)iD-P-3) , 



iD-P-2){D-P-3) dab 

2 ^D-P-l' ^^-^^1 

for any a,b = 1..P. One can then easily compute the conserved mass and the diagonal 
components of the tension tensor which are given by the relations ( p.29D below. Note that 
in the general case, the tension tensor has also the non-diagonal components: 

Recall now that we are interested in discussing the thermodynamic properties of black 
objects whose metrics are asymptotic to the form given in ( |2.9| ), which is appropriate when 
considering nonuniform extended objects with (p + l)-compact spatial directions. In order 
to accommodate this demand one simply has to replace in the above formulae P with p + 1 
and single out one of the coordinates, say the {p + l)-coordinate and denote it by z, while 
the p remaining compact directions will be denoted by Y^. Further note that for the ansatz 
considered in ( p. 91) one has Ca = Cy, dab = for a, 5 = l..p and Cp+i = c^, while Vp = VpL, 
where L is the length at infinity along the ^;-direction. 

3.2 Another counterterm choice 

For comparison and also as a check of the correctness of the results obtained in the previous 
section, we shall complement these results with those obtained by using the following 
boundary counterterm: 

= ^ / d^-'xV^Vn, (3.26) 
oT^*^ JdM 



where C(£) p) = ^J {D — p — ?,) / {D — p — 4)^ and TZ is the Ricci scalar of the boundary 
geometry. 

Varying the total action with respect to the boundary metric hij, we compute the 
boundary stress-tensor: 



where Kij is the extrinsic curvature of the boundary and we denote ^ = C(£) p) yfXjTZ. 

One should remark here that the counterterm choice is not unique, other choices being 



possible as well depending on the boundary topology (see |14, 25, 27, 20 1 for a related 
discussion). Our choice of using ( 3.26| ) was motivated by the fact that the general expression 
for the boundary stress-tensor is very simple^. 



^For example, we have verified that for D = 5, 6, 7, the same expression for the mass and tension are 



found by using a counterterm on the form (3.26), with replaced by Ti?^'^ j \^7i? —J^- However, in that 
case, the expression of the boundary stress-tensor is much more complicated than (3.27), see e.g. ref. p^ . 
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Using the counterterm ( 3.26| ) for geometries whose asymptotic form is given by (p.S|), 
to leading order term the relevant components of the boundary stress-tensor are given by: 

1 {D -p- 3)q -Cz- pcy 



Tt 



rpY, ^ 1 Ct-c,-{D- A)cy 

rr.Z _ 1 Ct-{D-P- 3)C^ - pCy 

Note that we already made use of the relations Ca = Cy, for a = l..p and Cp+i = c^. The 
mass and the gravitational tensions are the conserved charges ( |3.2|) associated to d/dt, 
d/dy^ and d/dz, respectively (note that d/dz is an asymptotic Killing symmetry of the 
metric). For a straight torus geometry along the brane directions it becomes apparent that 
the tension-tensor T"'' is diagonal. Keeping this in mind, when defining the gravitational 
tension T"" along one direction it is also more convenient at this point to extend the 
integration along all the spatial directions excluding only the direction y". For a given 
direction this amounts to multiplying T"" with the factor Vp/La- 

3.3 The global charges 

For both the counterterm choices discussed above, one finds the following expressions for 
the gravitational charges: 

M = ^"^^g^^^"' {{D-P- 3)Q - C. - pCy) , 

% = {ct-{D-p- 3)c. - pCy) , (3.29) 

where r2£)_p_3 is the area of the unit {D — p — 3)-dimensional sphere^. 
Usually in the literature one also defines a relative tension 

n/ = ^ , (3.30) 

(with / = z,i) which measures how large the tension along the /-direction (or binding 
energy) is relative to the mass. 

The solutions we construct will also have an electric charge which is computed accord- 
ing to 



Qe 



I ^.U.^Y^^O (3.3.) 

(with Qe defined in ( |2.11| )), and a dilaton charge which is read from the asymptotic 
expansion ( 2.12| ) of the scalar field 

Qd = ^^^^^(I)-p-%,. (3.32) 



"Note that these results are consistent with those previously derived (by other means) in the literature 
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As we shall discuss in the next Section, is not an independent charge for the solutions 
generated by using the Harrison transformation and it can be expressed in terms of mass, 
z-tension and electric charge. The same holds for the gravitational tensions along the 
y*-directions. 



4. KK static, vacuum black objects 'u\ d = D — p dimensions 

The charged black objects discussed in the next sections are constructed using as seed 
solutions the known vacuum static BHs and black strings va. d = D —p dimensions {i.e. with 
only one compact direction). In absence of matter fields (or of a cosmological constant), 
these configurations are trivially uplifted to D— dimensions by adding p— flat directions Y^. 
Therefore it is of interest here to briefly review their basic properties (see |28| for a more 
detailed overview of the KK solutions). 

Already in this case, the vacuum Einstein equations admit a variety of interesting 
solutions, approaching asymptotically the Ai'^^^ x background. Of interest here are 
the black strings and the black holes, which are distinguished according to their horizon 
topology. Both types of solutions are usually described within the following ansatz with 
three unknown functions^ A,B,C: 

dsl = -e'^^^'^^mdt' + e2^(^'^) (-^ + dz^) + e'^(^'^^R\r, z)dcol,, (4.1) 

where A,B,C vanish at large coordinate r. f{r) and R{r,z) are suitable background 
functions. Usually one takes /(r) = 1 — {ro/r)'^~'^, R{r,z) = r for black strings, in which 
case To > corresponds to the horizon radius. The usual choice for black holes is /(r) = 1, 
while R{r,z) has a complicated form (with R{r,z) ^ r in the asymptotic region), see 
Appendix A. 

Unfortunately the only solutions known in closed form (except for the bubbles, which 
we do not consider in this work) have no z— dependence and describe uniform black strings 
(UBSs) with A = B = C = (i.e. the direct product of the Schwarzschild-Tangherlini 
solution in (d — l)-dimensions with a spacelike circle). The existence of both caged BHs 
and NUBSs has been asserted so far mainly based on numerical results (note, however, 
that analytic approximations of the small BH solutions have been presented in |3^). Such 
solutions have been found at the nonperturbative level by solving a set of three partial 
differential equations for the metric functions A, B, C, subject to a suitable set of boundary 
conditions. The technical aspects involved in the numerical construction of these solutions 
are presented in Appendix A. 

The horizon topology of the black string solutions is S"^^^ x S^. As shown by Gregory 



and Laflamme |31|, a uniform solution is unstable below a critical value of the mass. As a 



result, there also exist nonuniform strings, which emerge smoothly from the UBS branch 



at a critical point, where its stability changes |31]. Keeping fixed the horizon coordinate 



'^An ansatz involving only two undetermined functions has been proposed in j2^], for a special choice of 
coordinates. However, it remains an interesting open problem whether one can use that ansatz for numerical 
calculations. 
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ro and the asymptotic length L of the compact direction, these solutions depend on a 
single parameter, specified via the boundary conditions and physically corresponding to 
the Hawking temperature of the solutions [^2|. Varying this parameter, the nonuniform 
strings become increasingly deformed (see Appendix A for details). A measure of the 
deformation of the geometry of these solutions is given by the nonuniformity parameter A 

m 

where T^max and T^min represent the maximum radius of a (d — 3)-sphere on the horizon 
and the minimum radius, being the radius of the 'waist'. Concerning the thermodynamical 
stability of the NUBS solutions, it turns out that, at least for d = 5, 6, the entropy decreases 
with temperature, i.e. they are thermally unstable. 

Apart from the black strings, KK theory also admits as solutions a branch of caged 
black holes with an event horizon of topology S'^~'^ . For the metric ansatz ( [4.1D , the 
event horizon is still located at a constant value of the radial coordinate, r = ro, with 
— oo < r < 0. (Note that the functions A, B, C contain in this case suitable 'background' 



parts that enforce the required behaviour on the boundaries and the topology difference 
between the event horizon and infinity, see Appendix A.) 

For small values of the mass, these solutions look like the Schwarzschild-Tangerlini 
solution in d dimensions. Once one increases their horizon radius for a fixed value of the 
asymptotic length L of the extra-direction, they deform as they start to 'feel' the compact 
direction and their mass increases. Analogous to the NUBS case, the geometrical prop- 
erties of the black holes can be described in terms of two parameters denoted here by Req, 
the equatorial radius of the horizon, and Laxis, the proper length along the exposed sym- 
metry axis. Again, the entropy of the caged BHs is a decreasing function of temperature, 
i.e. the considered BHs are thermally unstable^. 

Although a number of details remain to be clarified, the numerical results suggest that 
the black hole and the NUBS branches merge at a topology changing transition Q, p5| . 
This scenario is supported by all numerical data available. For NUBS, the normalized 
thermodynamic quantities change between one (the uniform solution) and some critical 
values. The normalized mass, relative tension and entropy of the caged black holes vary 
between zero and some critical values that seem to correspond to those where the NUBS 
branch stops to exist. At the same time, the Hawking temperature of the BHs varies 
between infinity (for vacuum state with M = T = 0) and a minimal value. The set of 
black hole and the set of NUBS solutions then appear to merge at critical values of these 
quantities. In a (n, M) or (n, T) diagram, this feature can be seen looking at the vacuum 
data {U = 0) in Figure 5. (For convenience, the thermodynamic quantities there are 
normalized by their value for the critical UBS solution.) This scenario also implies that 
the 'waist' of the NUBSs would shrink to nothing, locally forming a cone geometry. At 
that point, both Laxis and Rmin go to zero, while Req tends to Rmax- 



*Note however, that the recent results in [g4[ indicate the existence of a more compUcated picture, with 
the possible existence of caged BHs possessing a positive specific heat. 
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For both black holes and black strings, the following Smarr relation holds 

{d - 3)M - LT 



d-2 



ThS, (4.3) 



where Th is the Hawking temperature and S is the entropy of the solutions (which is one 
quarter of the horizon area). The first law of thermodynamics, which for the KK black 
objects takes the form: 

dM = TudS + TdL (4.4) 

is also satisfied. 

5. Charged black objects in KK theory 

The vacuum configurations discussed above can be uplifted to D-dimensions and used to 
generate solutions of the Einstein-dilaton-antisymmetric form theory by applying a Harri- 
son transformation, as described in Section 2. The asymptotic structure of the metric and 
the topology of the event horizon remain unchanged in the process. 

5.1 The solutions and the geometry 

In our approach, the two metric parts which enter ( |2.4D are 

gf{x)dy'dy^ = -e2^«(^'^)/(r)dt2 + j^idY'f , (5.1) 

i=l 

ap j(r) ^ 

where ^o(^5 z),BQ{r, z) and Co(r, z) are the functions of a vacuum seed solution. 

The simplest charged solution one can construct by employing the Harrison trans- 
formation described in Section 2, uses a vacuum uniform black string solution as seed 
configuration. The resulting solution corresponds to a (p-|- 1)— brane solution with the line 
element: 

ds' = - (i + (^)^-^-^T^)"'"^''"'''^ [ - (1 - rf)^-^-')dt' + ELMY'?] 

+ (1 + (^)^-^-^T^)2'^(^+^) + dz^ + r^du^h.,.,) , (5.2) 

and matter fields 



^ = _2a{D - 2)alog(l + (m)^-P-4^), 

with a given by (^^) and < U < 1. For any value of a, U, the surface r = tq is still an 
event horizon, while r = is a curvature singularity. As we have previously mentioned, 
the geometric structure outside the event horizon and also the asymptotic structure are 
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preserved by the solution generating procedure. The inclusion of the matter fields affects 
the geometry in both z— and y*— directions, however in different ways. One can see that 
the proper length of the z circle increases for any finite r, while the opposite holds for a 
y*— direction. Moreover, for a given value of r, the area r2£)_p_3 of the D — p — 3-sphere 
increases with U . 

The extremal limit of (^]^) is found by taking [/ — > 1 together with a rescaling of tq 
such that Tq = (l — f/2)c^~P-4 ^^^^ reads: 

ds^ = - (l + (^)^-P-4)-^"(^--^) l^-dt^ + ±{dYr^ (5.3) 



Aty,...y, = ^/2{D - 2)a (l + {^)D-P~A'' , = -2a(D-2)alog(l + (-)^-P-4), (5.4) 



c being a positive constant. 

The basic features of the nonextremal uniform configurations are shared by the charged 
solutions generated starting from NUBS and caged black holes. In both cases the geometry 
is described by the metric form: 

ds^ = _e2A(r,.)^(^)^^2 ^ ^2Bir,z) + dzA (5.5) 

\f{n J 

i=l 

where 

A{r, z) = Ao(r, z) - a{D - p - 3) log Q{r, z), B{r, z) = BQ{r, z) + a{p + 1) log r2(r, z), 
C{r, z) = Co(r, z) + a{p + 1) log n{r, z), Ci(r, z) = -a{D - p - 3) log n{r, z), (5.6) 

and 

^) = 13772(1 - C/'e2^°(^'^)/(r)). (5.7) 
The only nonvanishing component of the {p + l)-form and the dilaton are given by 

Aty,..y, = V2(D-2)ae2^o(^'^)/(r) c/> = -2a(I) - 2)a log J7(r, z). (5.8) 

Again, the position of the horizon is not affected by the Harrison transformation. Once one 
knows the expressions of the functions AQ{r, z), BQ{r, z), Co(r, z) for the vacuum solution, 
one can derive all the geometric properties of the corresponding charged configurations. 
However, in absence of an explicit analytical solution, it is unclear how to derive the 
extremal version of the black objects with an explicit z-dependence starting from the 
non-extremal solution (^.51), (^.^). Therefore, for the rest of this work we shall consider 
nonextremal solutions only. 
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5.2 The charges of the solutions 

A straightforward computation based on the generic asymptotic form of the metric func- 
tions ( 2.1C| ) leads to the simple relations: 

C/2 



ct = cf\l + 2a{D -p-2,) 



c, = 40)+2cl°^a(p + l) 



(5.9) 



cf\ ck^^ being the constants in the asymptotic expansion of the seed solution with U = 0. 
The expressions for the electric charge Qe and the dilaton charge Qd, as derived from 



(|2l0|),(|3:3l|),( p:3^) , are: 



U 



Qd 



levrG 



,(0) 

t ! 



- p - 4) 



SO) 



(5.10) 
(5.11) 



There is also a one-to-one correspondence between the solution generating parameter U 
and the electrostatic potential difference $ between the event horizon and infinity 



(5.12) 



As argued below, in a thermodynamical description of the system, $ corresponds to the 
chemical potential. 

One can see that the mass and tensions of a charged solution can be expressed as linear 
combinations of the corresponding quantities of the vacuum configuration: 



' m' 




Tyi 









Oil ai2 ai3 

021 022 023 
OSl ^32 033 



Mo 



(5.13) 



where 



On = 1 + 



021 



2{D - 2){D - p - 'i)a 

{D-p-2) l-U^' 
2{D - 2){D - p - 3)a 



013 



{D-p-2)Li 



023 



2{D - 2)aL [/2 
" ~(L>-iJ-3) l-[/2' 
2{D -2){D - p-3)aLi 

{D-p-2)L l-C/2' 



while ai2 = 031 = 032 = 0, 022 = 033 = 1. (Note that the tension along the z— direction 
does not change when adding an electric charge to the system). 

The electric charge can also be expressed as a function of mass and z— tension of the 
vacuum solution: 



with a similar relation holding for the dilaton charge. 



(5.14) 
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The initial vacuum configuration has two global charges, Mq and T^. The Harrison 
transformation introduces via the parameter U a new charge. This implies that Qd and 
% are not new global independent charges and therefore they do not enter the thermody- 
namical description of the system since they cannot be varied at will. Indeed, one finds 
the following quadratic relation 

Ql =fl- ^(^-5^-r» lg,. (5.15) 



Ki{M{D-p-3)-L%) + Qd V D-p 
with 

^ {D - p- 4) (D - p- 2 - 2{D - 2){D - p- 3)a)' ^ ' ' 

which determines the dilaton charge in terms of the mass, tension along the z— direction 
and the electric charge. The tensions along the y*-directions can also be written as 

D-p-3\ 2{D-2)U^a zj k ) 

5.3 The thermodynamics 

The Hawking temperature Th of a black object is computed either by requiring the absence 
of conical singularities on the Euclidean section (which is found by analytical continuation 
t — it) or from the surface gravity of the solutions in the Lorentzian section. As usual 
with black objects in Einstein gravity, the KK solutions considered in this work have an 
entropy which is one quarter of the event horizon area. 

One can easily prove that the relation between Hawking temperature Th and the 
entropy S of the charged solutions and the corresponding quantities ^//(o) and Sq of the 
vacuum seed solution is 

Th = Th^o){1-UY^''-^\ 5 = 5o(1-C/2)-"(^-2), (5.18) 

such that the product between the entropy and temperature remains invariant under the 
Harrison transformation. Thus adding an electric charge to a system with a fixed event 
horizon radius would decrease its temperature and increase the entropy, which agrees with 
the experience based on Reissner-Nordstrom solution. 

Considering the thermodynamics of these solutions, they should satisfy the first law of 
thermodynamics 

dM = THdS + ^dQe + %dL. (5.19) 

One may regard the parameters S, Qe and L as a complete set of extensive parameters 
for the mass M(S,Qe, L) and define the intensive parameters conjugate to them. These 
quantities are the temperature, the chemical potential and the z-tension. 
For any type of charged black object, the following Smarr relation holds 

D — p — 2 D — p — 2 
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which results from (|4.3| ) together with (|5.13| ), ( |5.14 ) and ( 5.18 ). 

In the canonical ensemble, we study solutions holding the temperature Th^ the electric 
charge Qe and the length L of the extra-dimension fixed. The associated thermodynamic 
potential is the Helmholz free energy 



F[Th, Qe, L]=M- ThS, dF = -SdTu + ^dQe + T,dL . 



(5.21) 



The situation of KK black objects in a grand canonical ensemble is also of interest, in 
which case we keep the temperature, the chemical potential and the tension fixed. In this 
case the thermodynamics is obtained from the Gibbs potential 



G[Th,<^,%] =E-ThS- $Qe - %L. 



(5.22) 



By using the relations ( 5.13| ), ( |5.14 ), one finds that both F and G can also be expressed in 
terms of the global charges of the vacuum seed solutions and the generating parameter U . 

It is well-known that different thermodynamic ensembles are not exactly equivalent 
and may not lead to the same conclusions regarding the thermodynamic stability as they 
correspond to different physical situations (see e.g. |3(:], [37|). Mathematically, thermody- 
namic stability is equated with the subadditivity of the entropy function. In the canonical 
ensemble, the subadditivity of the entropy is exactly equivalent to the positivity of the 
specific heat at constant {Qe,L) 



Cr 



Th 



dS 



> 0. 



(5.23) 



Remarkably, one can prove that the specific heat of a charged configuration can be expressed 
in terms of the specific heat and entropy of the vacuum seed solution^: 



Cq,l 



1 



Cq l(1 + - MD - 2)C/2) - 2a{D - 2)U'^Sq 



(1 _ f/2)2a(D-2) 



C/2 + i + 2a(Z)-2)C/2(l + 



Q,L - 



(5.24) 



So 



The condition for thermodynamic stability in the grand canonical ensemble is [37|, |38]: 

Cq,l > 0, eT ^ (II) > 0, (5.25) 

where e-r is the isothermal electric permittivity, which measures the stability of solutions 
to electrical fluctuations. This quantity can also be expressed in terms of the quantities of 
the vacuum seed solution^" 



+ l + 2a{D -2)U^{1 + 



Sq 



(5.26) 



relation (ksl) implies also cf^ ~ T„^o)So- Thus, one finds rf/(o)S'od((7/(l-l7'')) + f//(l-;7'')d(rH(o)5'o) = 



In evaluating Cq,z, one uses (|5.1q) together with the expression (5.10) for the electric charge. The Smarr 



and substitutes back in (5.23). 

^°Here one uses the relation = dTH(o) - '2a{D - 2)U/{1 - U'^)TH{o)dU , which follows from ( ^ 
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We have seen in the previous Section that the specific heat of the vacuum solutions 
Cq 2, is a strictly negative quantity for both black holes and black strings. However, one 
can see from (|5.24 ), ( p.26 ) that Cq^l and e-r may become positive for some range of U, 



a, p. 

Here we would like to point out that for any black object constructed by using the 
d = 5,6 vacuum solutions discussed in Section 4, the condition ( p.25| ) for thermodynamic 
stability in the grand canonical ensemble cannot be satisfied. To prove that, we express 
Cq^l and er as 

Cq,l = (,_^.)L(o-2) .go (i+fJ^TTy - l) ' = + U^){1 + «(x + 1)),(5.27) 

where we denote x = Cq i^/Sq, a = a {D - 2)2C/V(1 + U^)- Moreover, we have verified 
that Cq ^ + Sq < for all vacuum seed solutions with d = 5, 6 {e.g. this quantity is 
Cq ^/ Sq + I = —{D — p — A) for UBSs). Therefore x + 1 < and, as a result, ex and Cq^l 
cannot be both positive at the same time. (This is a general result, valid for any values of 
a,p, U .) Thus we conclude that all charged KK solutions generated from d = 5,6 vacuum 
solutions are unstable in a grand canonical ensemble. 

Of course, this does not exclude the existence of stable solutions in a canonical ensem- 
ble. In the next subsection we shall find such configurations for both black holes and black 
strings (i.e. with p = 0). A similar analysis for a generic p does not seem to be possible, 
except for the case of uniform brane solutions (thus without an explicit dependence on 
z), which we shall consider next. The quantities which enter the thermodynamics can be 
derived^^ easily in this case from (^.2|) 



M = ^""-^^^f^ {{D-P- 3)(1 - U') + 2{D-2){D-p- A)U'a) f-^, (5.28) 

_ nD-p-4Vp D-p-4 ^ _ np^p^^VpL /— — U D-p-4 

levrG ' levrG V^{D - ^KD - p - 4)^—^r, 

where ro is the event horizon radius. One can see that the charge to mass ratio of these 
solutions has a nontrivial dependence on a,p,U: 

Qe ^ ^2iD-2)a{D-p-A)U 

M {D -p-^){l-U^) + 2{D -2){D -p-4)U^a ^ ' ' 

The specific heat and isothermal electric permittivity of the uniform solutions are given by 

2a{D -2){D-p- A)U^ - (l + [/2)(Zj - p - 3) 
^Q,L (1 _ f;2)2a(D-2) I + u2 _ 2a{D - 2){D - p - A)U^ ' ^ ^ ^ 

= + 1 - MD -2){D-p- 4)C/2) . (5.31) 

For given a,p this implies the existence of a phase transition for a critical value of U 

Uc = ^ (5.32) 

^2a{D - 2){D - p - i) - 1 



D-p-i 



^'^For the solutions with a dependence, we have used the numerical data in [Q, [ p9[ . 

^^The thermodynamics of uniform black strings and p— branes has been discussed in (see also jjl] ] 



and references therein). 
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where Cq,l changes the sign. All solutions with U > Uc (thus with a large enough electric 
charge) are stable in a canonical ensemble. Also, one can easily show that Cq^l is always 
negative for large enough values of a (this property is shared by the solutions with a 
dependence on the extra-dimension z). 

This issue is especially interesting in connection with the Gubser-Mitra (GM) conjec- 



ture 1 42], that correlates the dynamical and thermodynamical stability for systems with 
translational symmetry and infinite extent. According to this conjecture, the appearance 
of a negative specific heat of a black brane is related to the onset of a classical instability. 
Supposing that this conjecture holds also for the charged solutions discussed in this paper, 
one would expect all charged, nonextremal black branes to be unstable, irrespective of the 
value of the electric charge^^. This is indeed the case when considering charged solutions 
in a grand canonical ensemble. This agrees also with the results based on the map ( 5.13| ). 



As shown in |31], the vacuum uniform solutions with a given L are unstable below a critical 
mass A consequence of the Harrison transformation is that the unstable mode of the 

vacuum solutions would be present also for the nonextremal charged solutions. Thus, for 
a given U, one expects a branch of charged, nonuniform solutions to emerge for a critical 

2{D~2){D~p-i) - 
(D-p-3)(l-C/2) I 



value of the mass of the uniform solution Mc = Mc'^\l + ^rrtl U^a) 



5.4 The p = case: charged black strings and caged black holes 

To simplify the general picture, we shall analyse for the rest of this section the simplest 
case of charged nonuniform black strings and black holes in D = 5, 6 KK theory, i.e. 
with a single extra-direction, p = 0. As we shall see, the solutions have a rich structure 
already in this case. 

Starting with the geometrical properties, we note that for the approach in this work. 



the charged BHs and the NUBSs are described by the ansatz (4.1) (we recall that /(r) = 
1 ~ {i"o/f)'^~^, R{^, z) = r for black strings, while /(r) = 1 and i?(r, z) has a complicated 
expression for caged BHs, see Appendix A). As seen from ( ^.61 ), the inclusion of an electric 
charge changes the shape of the functions A,B,C. 

In Figure 1 we exhibit the metric functions —gtt/f = e^"^, Qzz = e^"^, and gge/r"^ = e^*^ 
for a typical D = 6 vacuum NUBS solution (U = 0) and a charged one with U = 0.85 
in Einstein-Maxwell theory (a = 0). The nonuniformity parameter of these solutions as 
defined by (|4|) is A = 0.65. 

To obtain a more quantitative picture of the NUBS metric functions, we exhibit A, B, 
and C in Figure 2 for three values of U and two fixed values of z. The solutions there have 
a dilaton coupling constant a = 1 and A = 0.93. The corresponding matter functions for 
the same configurations are shown in Figure 3. One can see that although the shape of the 
functions is qualitatively the same in both cases, the precise values taken by A, B, C in the 
region close to the horizon are very different. However, in both cases the metric functions 
exhibit an extremum at the horizon for z = L/2. 



^^The Gregory-Laflamme instability of D = 5 electrically charged black strings has been discussed in the 
recent work [ p3| . 

^"^These are the only cases where data is available for both the BHs and NUBSs branches. 
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Figure 1. The metric functions , , and e^^ of a typical D = Q nonuniform string solu- 
tion arc shown as functions of the coordinates x — ^Jr"^ ~ j ^ f'^'' ^ vacuum solution (left) and a 
charged solution with U = 0.85 (right). (Note the different scale of the fmictions in these cases.) 

The proper length of the z— circle increases with U, for any finite value of U. (The in- 
put parameters of the solutions in Figures 1-3 are ro = 1, L = 4.9516). 

An interesting point here concerns the deformation of the horizon of the NUBS solu- 
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tions. One can see from ( |5.5D that although the charged solutions have different values of 
Rmax, Rmin from the vacuum case, the nonuniformity parameter A as defined in ( |4.2| ) does 
not depend on the value of the chemical potential (or of the electric charge). 
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Figure 2. The metric functions A, B, and C oi D = 6 charged nonuniform black string solutions 
are shown as functions of the radial coordinate x = y'r^ — r^, for three values of the parameter U. 
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Therefore, the quahtative features of the horizon geometry remain the same as in the 
vacuum case. In particular, the parameter A approaches very large values for solutions 
constructed using seed NUBS close to the transition point (this occurs for any value of U). 



z/L=0.23 



0.8 



-0.4 





A, 






■ 








U=0 






U=0.964 






U=0.761 









0.5 0.75 



0.8 ^ 

0.4 - 

— 

-0.4 - 

-0.8 -/ 

-1.2 - 

-1.6 - 

-2 - . 

-24 \[_ 




z/L=0.46 



u=o 

U=0.964 
U=0.761 



Figure 3. The dilaton <j> and the electric potential At are shown for the same solutions exhib- 
ited in Figure 2. 



The geometrical properties of the caged BH solutions can be discussed in a similar 
way. In Figure 4 we exhibit the metric functions e^^, e~^^ and e^*^ for a typical D = 6 
vacuum BH solution (U = 0) and a charged one with U = 0.85 in EM theory (a = 0). 
The position of the event horizon there was = — 1 (see Appendix A for more details 
on the seed solutions). As a generic feature, for any value of the parameter U, the ratio 
between Req and Laxis remains the same as in the vacuum case. At a technical level, this is 
a consequence of the fact that the Harrison transformation ( p.5[ ) involves just a conformal 
transformation of the spacelike part gapdx'^dx^ of the seed geometry. 

This observation is relevant in connection with the issue of a topology changing transi- 
tion between the black hole and the NUBS branches. Since, from ( |5.13| ), adding an electric 
charge implies just a recombination of the global charges of the system, we argue that the 
picture one finds here would be qualitatively similar to that valid for vacuum solutions. 

To clarify this point, we show in Figure 5 the mass M and temperature T versus the 
relative string tension n, for the nonuniform string branch as well as for the black hole 
branch^^, for several values of the parameter U. The black hole data are taken from p5| . 
(For convenience, we normalize in this subsection the thermodynamic quantities by their 
values for the corresponding charged UBS with the same value of U.) 



^^The same picture is found when considering the (n, M) diagram in D = 5 and the (n, T) plot in D = 6, 
respectively. For all values of U, we note qualitative agreement of the shape and the relative position of 
the charged nonuniform string branch and the charged black hole branch in five dimensions with the shape 
and relative position of the corresponding branches in six dimensions. The backbending noticed in for 
the vacuum case is present also for charged NUBS, both in five and in six dimensions. 
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Figure 4. The metric functions e^"^, e , and e^*^ of a typical = 6 KK black hole solution are 
shown as functions of the coordinates r, z for a vacuum solution (left) and a charged solution with 
U = 0.964 (right). 

One can see that we have all reasons to expect that the topology changing transition 
scenario proposed in the vacuum case, holds here also. However, as seen in Figure 5, a 
nonzero electric charge translates into a different position in the (M, n) plane (and {Th, n), 
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Figure 5. Left: The mass M of the D — 6 charged nonuniform string and black hole branches 
is shown versus the relative string tension n. Right: The Hawking temperature Th of the D — b 
charged nonuniform string and black hole branches is shown versus the relative string tension n 
(these plots are given for three values oiU). M , n and Th correspond to the quantities normalized 
by the corresponding values of the uniform string solutions with the same values of U. 

respectively) of the putative transition point. We have also noticed that as U increases 
towards its maximal value, the M(n) curves converge to a critical one, which is close to 
that corresponding to U = 0.964 in Figure 5. The same holds for the (n, Th) diagram. 

The general thermodynamical properties of the charged solutions are rather compli- 
cated as compared to the vacuum case, with a nontrivial dependence on the dilaton coupling 
constant a. However, the qualitative features are again the same for both D = 5 and D = 6 
solutions. 




Figure 6. The scaled isothermal electric permitivity is shown as a function of the normalized 
temperature and the parameter U for the D = 6 charged caged black holes (left) and nonuniform 
black string strings (right). 
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Figure 7. The scaled isothermal electric permitivity is shown as a function of the dilaton coupling 
constant a for the D — 5 charged caged black holes (left) and nonuniform black strings (right). 
The plots here correspond to a fixed normalized temperature Th — 0.7 and several values of the 
chemical potential U. 



Also, as argued above, both the BHs and the NUBSs are thermodynamically unstable 
in a grand canonical ensemble^^. Here we start by shovi^ing in Figure 6 the scaled isothermal 
electric permitivity-'^'^ 



So ' 



(5.33) 



for caged BHs and NUBS solutions as a function of the normalized temperature and the 
parameter U. These plots correspond to a value a = of the dilaton coupling constant, 
i.e. the EM theory. One can see that, for NUBS, €t is a decreasing function of U for any 
allowed value of the temperature. However, it becomes negative for large enough values of 
the chemical potential. The NUBS solutions close to the transition point have a positive 
€t for small U only. As seen in Figure 6 (left), the picture is more complicated for caged 
BHs. The small charged black holes {i.e. with large values of Th) have always a positive 
€t- However, larger black holes have a negative isothermal electric permitivity, for large 
enough U. 

This picture depends however on the value of the dilaton coupling constant a. Our 
study indicates that for both types of black objects ex becomes positive for large enough 
values of a, irrespective of the values of U, Th (see Figure 7). 



^® As discussed in , an unexpected feature of the NUBS solutions is the occurence of a new branch at 
a critical (minimal) value of the relative string tension, close to transition point. Along this branch (which 
has a small extension in the (M,n)— plane, see Figure 5), the mass decreases with n, while the specific heat 
is still negative. In what follows, to simplify the general picture, we shall ignore this secondary branch for 
all plots of the specific heat and isothermal electric permitivity. The same holds for caged BH solutions, the 
last points for the data in |^] indicating the existence of a (small) secondary branch close to the transition 
point, see Figure 5 (this feature was confirmed by the recent results in js^]). This secondary BH branch 
was also ignored when discussing Cq,l and er- 

^"^It has the same sign as the isothermal electric permitivity and stays finite as t/ — > 1. 
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Figure 8. The specific lieat is sliown as a function of U for tlic D — 5 charged caged black 
holes (left) and nonuniform black string strings (right). The plots here correspond to a vanishing 
dilaton coupling constant a — and several values of the normalized temperature. 
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Figure 9. The specific heat is shown as a function of the scaled temperature for the D ^ 6 
charged caged black holes (left) and nonuniform black strings (right). These plots are for a fixed 
value of the chemical potential C/ = 0.5 and several values of the dilaton coupling constant. 



This can be understood by observing that for large o, the term becomes dominant 

in the numerator of the isothermal electric permitivity. 

In Figure 8 we present the specific heat as a function of the parameter U for several 
values of the rescaled temperature for black strings and caged black holes in D = 5 EM 
theory (a similar picture is found in D = 6). The specific heat of the black strings is 
always negative for small values of U. However, for any Th, it changes the behaviour for a 
critical value of the chemical potential U = Uc (which is fixed by the condition eT{Uc) = 0). 
Cq^l exhibits a discontinuity at Uc, indicating the existence of a phase transition in the 
thermodynamical system. The value of Uc decreases as the scaled temperature increases. 



-28- 



BH 



NUBS 





1 U=0.9 




ii ; ; i\ u=o.7 




M / i U=0.6 




/ i ! \ / : U=o.5___ 




i u=o.4 









450 



O 



. U=P,9 .. 
. .U.=Q.7 . 

.11=0.5 

. U=P..6 .. 
...U=M-. 

U=0.2 



Figure 10. The specific lieat is shown as a function of the dilaton couphng constant for the 
D — 6 charged caged black holes (left) and nonuniform black strings (right). The plots here cor- 
respond to a fixed value of the scaled temperature Th — 0.8 and several values of the chemical 
potential. 

The picture is slightly different for caged black holes. The scenario we found for NUBSs 
applies only for solutions close enough to the critical region where the topology changing 
transition is expected to occur. The solutions generated from small caged BHs with a large 
enough scaled temperature have a negative specific heat (see Figure 8 left). 

Again this picture depends on the value of the dilaton coupling constant. In Figure 9, 
the specific heat is shown as a function of the scaled temperature for D = 6 solutions with 
a fixed value of U and a range of a. 

The dependence of the solutions on the dilaton coupling constant a is further illustrated 
in Figure 10, where we show the specific heat as a function of a for several values of U. 
For small values of a, one notices the possible existence of configurations with a positive 
specific heat, for both caged BHs and NUBSs. One can see that for large enough values of 
the dilaton coupling constant, the specific heat becomes negative irrespective of the value 
of U, for both black holes and black strings. This can by understood by noticing that, 
since a — > for large enough a, the sign of the specific heat ( 5.24| ) of the charged solutions 
is fixed by Cq j^, which is a negative quantity. 



6. Slowly rotating charged NUBS 

It would be interesting to generalize the charged KK configurations discussed in the pre- 
vious sections by adding a nonvanishing angular momentum. Unfortunately, there is no 
simple procedure to spin a given static system. Moreover, all known techniques fail for the 
case of KK solutions with a dependence on the extra-dimension. Some of the reasons for 
this are that the great majority of the solution generating techniques rely heavily either 
on high degrees of symmetry in the metric (such as the geometries in the generalized Weyl 
class) or on the possibility of dimensional reduction down to 3-dimensions where target 
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space symmetries can be used to add rotation to a given static solution. As we have pre- 
viously stated, solutions with a dependence on the extra-dimension clearly do not belong 
to these classes of geometries and therefore any rotating non-uniform solution has to be 
found by brute force solving of the field equations. 

To the best of our knowledge, the only rotating solutions with a dependence on the 
extra-dimension correspond to the NUBSs discussed in [44|. The solutions there have been 



constructed in D = 6 dimensions, for the special case of two equal magnitude angular 
momenta. The angular dependence factorizes in this case, and the problem reduces to 
solving a set of partial differential equations with dependence only on the radial variable r 
and the extra-dimension variable z. D = 6 charged nonuniform black string configurations 



in heterotic string theory were also constructed in |44] by employing a solution-generating 
technique. It would be interesting to extend these results for other values of D and non- 
equal angular momenta. The case of a black string with a single angular momenta is of 
particular interest, given the complicated picture one finds already for D > 4 asymptotically 
Minkowski black hole solutions ||4^. However, this leads to a difficult numerical problem, 
since the equations would depend on three variables. Moreover, as far as we know, rotating 
EM-dilaton black strings do not exist in closed form even for the uniform case^^, and one 
has to rely on numerical methods to construct them. 

As an extension of the results in the previous sections, we present in what follows a 
numerical study of the NUBS solutions with a single angular momentum, in the limit of 
slow rotation (the case of UBSs is discussed in Appendix B, where a closed form solution 
is presented). For such configurations, the equations depend only on the radial variable r 
and the extra-dimension variable z and are solved to first order in the angular momentum 
parameter. To simplify the general picture, we shall restrict here to the case of EM theory, 
although the results can easily be generalized to include a dilaton field. 

In this approach, one starts with the following general ansatz depending on r, z and 
the angular variable 9, describing charged rotating NUBS solutions with a single angular 
momentum: 

+ e^^^*-'^'^) sin2 9id^ + J(r, z, 9)dtf + e^^^^'^'^) cos^ 9divl_ 

A = A^{r,z,9)dxf', (6.1) 

(we recall that /(r) = 1 — (ro/r)-^"'^ in this case). Then we suppose the following expansion 
of the functions in ( |6.1| ) in terms of a small parameter a (which should not be confused 
with the dilaton coupling constant): 

Xi{r, z, 9) = Xf\r, z) + aXf^ (r, z, 0), 

with Xf ^ (r, z) being the functions which enter the static solution {e.g. Air, z, 9) = A{r, z) + 
aAi{r,z,9), B{r,z,9) = B{r,z) + a Bi{r,z,9), C{r,z,9) = C{r,z) + a Ci{r,z,9) etc.). 

Except for special values of the dilaton coupling constant where the solutions can be generated from 
the Myers-Perry black holes following the techniques in |4q|, |47|]. 
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A straightforward but cumbersome computation shows that for infinitesimal a, the 
only term in the metric that changes is gttp and we find the following consistent metric 
form: 

\J[r) J 

- 2a sin^ ej{r, z)dipdt, (6.2) 

the angular dependence of the nondiagonal component being factorized. Similarly, the only 
component of the vector potential that changes is A^, the corresponding expression of the 
Maxwell field being 

A = At{r, z)dt + a sin^ Bh{r, z)d(p. (6.3) 

Thus the problem is reduced to solving the equations for j{r, z), h{r, z) in the background 
of the static charged solutions. These equations read (to simplify the relations, we denote 
here At = V): 

j,rr + ^ - {Vrhr + ^) " {A,rJ,r + ^) + 2j{A,rC,r + ^) 

+ {D- 5){Crj,r + ^) - 2j{Crr + ^) ' ^(D - 3)j{Cf, + ^) (6.4) 

+ ^(^.-2(Z)-3)C.-^) + i^^ = 0, 
r ' ' r r 

hrr + ^ + {A,rhr + ^) + {D - h){C,rhr + ^) + ^(2j(KrC. + ^) 

■ , V,zj,z,^ , 2e-^-^j\ , , ,^-5 , /\, 2(U-4)e^«-^g , 

- (T<r-J,r + — p)) + —^r^^^r + (^— + j)hr = . (6.5) 

The equation for j has a total derivative structure which helps us relate the total angular 
momentum to event horizon quantities: 

1 !^ d.e-^+(^-3)^r^-3 (^;, - f - 2jCr - hV,r) = lim,^oor^-3 (j,. - f ) .(6.6) 

For the slowly rotating solutions, the horizon radius is still determined by the equation 
/(ro) = 0. Also, the surface gravity and area of the event horizon do not change to order 
0{a). However, the angular momentum J and the magnetic dipole moment are evinced 
already at this order. Namely, the angular momentum is related to the asymptotic form 
of the metric by 

J sin^ 9 

^*^-"20z,-3Lr^-4' ^^-^^ 

while the magnetic dipole moment /x is read from the asymptotic form of the magnetic 
gauge potential h{r, z) 
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Figure 11. The functions j{r,z) and h{r, z) are shown for a typical D = 6 slowly rotating 
NUBS in Einstein-Maxwell theory (with x — \rr^—r^)- The parameters of this configuration 
are tq = 1, i = 4.9516, A = 0.47, U = 0.833. 

In practice, the idea is to start with a known vacuum NUBS solution and to charge 
it by using the Harrison transformation described in the previous sections. This generates 
static charged NUBS solutions of the EM equations {i.e. a = 0, no dilaton field). Then the 
eqs. ( |6.5| ) are solved with the background functions A, B,C. 

We have performed this computation for D = 6 NUBS solutions^^ with several values 
of the nonuniformity parameter A (with A defined by the equation (|4.2| )) and a range of 
U. The equations ( |6.4D were solved by using the same methods as those described in 
Appendix A for the vacuum NUBS case. In this scheme, one introduces a new radial 
coordinate x = y^r"^ — t'q such that the event horizon is located at x = 0. Then the 
equations ( |6.5| ) are solved subject to the following set of boundary conditions: 

j\x=o = jHe^'^^'^'^'K dxh\^=o = 0, j\x=oo = /iU'=oo = 0, dzj\z=o,L = 5^/iU=o,L = 0, (6.9) 

with ju a constant fixing the angular velocity of the horizon, Qh = ajn /tq . 

The profiles of the functions J, h for a typical solution are presented in Figure 11. One 
can see that they present a strong dependence on the extra-coordinate z, possessing a peak 
around z = L/2. 

A quantity of interest here is the gyromagnetic ratio g which is defined in the usual 
way by 

In D = 4 classical electrodynamics, the gyromagnetic ratio of a charged rotating body is 
always g = 1. However, the gyromagnetic ratio of an electron in Dirac theory is g = 2, which 



^^We expect slowly rotating NUBS solutions to exist for any value of D. However, the case I? = 6 is 
more convenient for numerics when dealing with KK solutions with a dependence on the extra-dimension, 
see e.g. the comments on that in [|48|, [p3[. 
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is precisely the value found for a rotating Kerr-Newmann black hole. When considering 
the higher dimensional version of these EM black holes^'^, the lowest order perturbation 
theory gives for the gyromagnetic ratio the result g = D — 2, which seems a natural 
higher dimensional generalization of the gyromagnetic ratio in D = A dimensions [^] , |^0[ . 
However, numerical calculations revealed that in higher dimensions, in a nonperturbative 
approach, the gyromagnetic ratio is not constant, but deviates from D — 2 for finite values 
of the charge |51|, |52], which was also confirmed by higher order perturbation theory p3|. 
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Figure 12. The gyromagnetic ratio g is shown as a function of U for several values of the nonuni- 
formity parameter A for D = 6 slowly rotating black string solutions in Einstein-Maxwell theory. 



However, so far most of the studies in the literature were limited to the case of black 
holes in asymptotically flat spacetime, the gyromagnetic ratio of solutions with compact 
extra-dimensions being not explored. As found in Appendix B, g varies for uniform black 
string solutions between D — 3 (no charge, U = 0) and ^^j^^^J"^^^^ (in the extremal limit, 
U ^ 1). One can see in Figure 12 that the picture is much more complicated for D = 6 
solutions with an explicit dependence on the extra-coordinate. For small values of A, the 
gyromagnetic ratio g still decreases with U. However, the situation changes for large enough 
values of the nonuniformity parameter A. Our numerical results suggest that, for any A, in 
the extremal limit, the gyromagnetic ratio of NUBS solutions approaches the same value 

gc = 32/11. 

It would be interesting to perform a similar computation for caged black holes and to 
consider the issue of topology change for spinning solutions. 

7. Conclusions 

The charged solutions in string theory and M-theory are usually obtained from the vacuum 
KK black objects in d-dimensions by sequences of boosts and dualities ||8|. This gave a 
precise connection between the rich phase structure of KK black objects and that of the 

^"Note that the D > 4 counterparts of the Kerr-Newmann black hole are not known in closed form. 
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corresponding branes in string/M-theory and provided predictions for the thermodynamics 
of supersymmetric Yang-Mihs theories on a circle. Our proposal here was to use instead 
the Harrison transformation derived in . This construction is more general since it holds 
for any value of the dilaton coupling constant a and any value of p (the brane extra- 
dimensions), in particular for the Einstein-Maxwell (EM) theory (p = 0). 

The main purpose of this paper was to use this technique in order to investigate the 
issue of charged black objects in KK theory for an arbitrary value of o, by using the known 
vacuum NUBS and the KK localized BHs as seed solutions. We have observed that both 
the geometrical and thermodynamical properties of the charged solutions can be derived 
from the corresponding D— dimensional vacuum configurations. Asymptotically flat black 
objects generated in this way from vacuum fields have a geometrical structure similar to 
that of the vacuum seed solutions. 

A natural question that arrises at this point is how to compute the various conserved 
quantities and the total action of the charged black objects. One particularly powerful 
approach to this problem is given by the counterterms method. For asymptotically flat 
spaces there have been proposed many such counterterms and their form relies heavily on 
the topology at infinity [p^]-|27|. However, in certain situations, these counterterms 
fail to regularize the gravitational action. One recent proposal put forward by Mann and 
Marolf seems to avoid these problems for asymptotically flat spaces. A secondary aim of 
this paper was to propose a direct extension of the Mann-Marolf counterterm to Hransverse 
asymptotically flaf spaces in KK theory. In Section 3 we showed that the total action and 
conserved charges can be expressed in a novel way in terms of the electric part of the Weyl 
tensor and that, moreover, the final results are consistent with those previously derived in 
the p-brane literature. 

The thermodynamics of charged black objects in KK theory has been considered in 
Section 5. There we have argued that these solutions are thermodynamically unstable in a 
grand canonical ensemble. Also, our results show that the thermodynamics of these black 
objects is affected by the value of the dilaton coupling constant. A general discussion of 
the charged black strings and black holes for the simplest case of a KK theory in D = 5, 6 
dimensions with a single compact extra-dimension has also been studied there. These 
solutions display a rich phase-structure already in this case, some of them being thermally 
stable in a canonical ensemble, for large enough values of the electric charge. 

The charged solutions we have discussed are likely to possess spinning generalization. 
Unfortunately, the Harrison transformation considered in our work is restricted to the case 
of static configurations. In Section 6 we have presented a numerical study of the uniform 
and NUBS solutions with a single angular momentum, in the limit of slow rotation. Our 
results provide evidence for the existence of such solutions and show that the gyromagnetic 
ratio depends on the nonuniformity parameter of the black strings. 

Although we have restricted our study to the (physically more relevant) region outside 
the horizon, it would be interesting to consider also the interior of these charged black 
objects. In the vacuum case, the ref. |54| presents numerical evidence, that the singularity 
touches the horizon as the horizon topology changing transition is reached. We expect that 
the inclusion of the electric charge would lead to a more complicated picture, at least in 
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the absence of a dilaton field. 

In this paper we have considered the simplest case of a single gauge field. In principle, 
one can easily add secondary gauge fields to the action (^]^) by considering sequences 
of boosts and dualities (this would also generate supplementary scalars, with a specific 
coupling to the gauge fields). The methods in this work can easily be extended to that 
case. On general grounds, we expect the features we have found to be generic; in particular, 
all solutions are likely to be unstable in a grand canonical ensemble. 

Also, as an avenue for future research, it would be interesting to study KK solutions 
with more general matter fields than the model The case of configurations with 

nonabelian gauge fields is of particular interest, given the rich structure of these systems and 
the interplay in this case between internal group symmetries and the spacetime symmetries 

m. 
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A. The construction of vacuum Kaluza-Klein solutions: technical aspects 
A.l NUBS solutions 

The NUBS solutions were constructed by several different groups [^, jS^, with 

different numerical methods, and are known for all dimensions up to eleven. 

The basic steps in the nonperturbative construction of these solutions are as follows. 
First, the numerics is performed by using a new radial coordinate x, which is defined such 
that the r, z part of the metric becomes conformally fiat^^, i.e. 

/djT 

and 

dsl = -e^^^''''^fir{x))df + e2^(^'^) {dx^ + dz^) + e^^^''^^\'^{x)dujl_^. (A.2) 

The constant xq in the above relation is fixed by imposing that the horizon is located at 
X = 0, where /(rg) = (moreover, x/r — > 1 in the asymptotic region). 

The Einstein equations = 0, G% + Gl = Q and = (where 9 denotes an angle 
of the D — ?) dimensional sphere and G\ the Einstein tensor) yield three elliptic equations 
for the metric functions A, B, C (these equations are given e.g. in Section 2 of ref. [0g|]). 
The remaining Einstein equations = 0, G% — Gf = yield two constraints. As first 
observed in I^Sj, these constraints satisfy Cauchy-Riemann relations 

a, {^gGD + {^g\{G% - GD) = 0, 9, (V^Gf) - {^g\{G% - GD) = 0.(A.3) 

^^For example, in D — x = ^/r^—r^. 
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Thus the weighted constraints satisfy Laplace equations, and they are fulfilled, when one 
of them is satisfied on the boundary and the other at a single point ||4^ . 

All known KK solutions with a dependence on the extra-dimension have the reflection 
symmetry w.r.t. z = L/2. Thus, in numerics it is enough to consider only half of the 
2;-domain. 

The NUBS solutions are then constructed subject to the following set of boundary 
conditions 

= C\ =0, 

D \X = 00 ' 

= 4, „ = a,C| „ = 0, (A.4) 



^1 

1 x= 


= B\ 


\x= 


=0 \x 




2=0,L/2 



lx=0 la;=0 

where the constant do > is related to the Hawking temperature of the solutions. Regu- 
larity further requires that the condition dxB\^_^ = holds for the solutions. 
The Hawking temperature and entropy of the NUBS solutions are given by 

Th = e^o-^»rO , S = Soy t e^^+^^-^'^^^'d^, (A.5) 

Jo 

where T^, Sq are the corresponding quantities of the uniform solution with A = B = C = 

T°=^, So = ^mn-srt', (A.6) 

and Aq{z), Bq{z), Co{z) are the values of the metric functions on the event horizon r = tq 
(x = 0). The mass and tension of these solutions are given by ( p. 29 ), where the constants 



ct, Cz are read from the asymptotics of the metric functions A and B, respectively. 

All known NUBS solutions were found by solving numerically the set of three coupled 
non-linear elliptic partial differential equations for A, B, C, subject to the above boundary 
conditions. For a better numerical accuracy, ref. proposed to employ a new radial 
coordinate f = x/{l + x), which maps spatial infinity to the finite value f = 1. 

The branch of nonuniform strings is obtained by starting at the critical point of the 
uniform string branch and varying the boundary parameter do, which enters the Eq. ( [A.^ ). 
The resulting solutions have a nontrivial dependence on z, the parameter A as defined 
by the equation ( |4.2| ) increasing with do for the main NUBS branch. Usually, the NUBS 
solutions are found with good accuracy (as verified e.g. by using the Smarr law (|4.3|)), 
except for the configurations close to the topology changing point (i.e. for very large values 
of the nonuniformity parameter A). 

A. 2 Caged black holes 

Given the mismatch between the horizon topology and the asymptotic structure of space- 
time^^ the numerics is much more difficult in this case, the only nonperturbative results 
available so far in the literature being for D = 5,6 dimensions [^], |61]. The ansatz 



^^It is interesting to notice the analogy with the D > A asymptotically Minkowski black holes with a 
nonspherical topology of the horizon which present the same feature. For example, the horizon topology 
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( |4.1| ) corresponds to that used by Kudoh and Wiseman in |3^, |35]. The independent 
results in |^l| were found by using a more complicated approach, with two different coordi- 
nate patches: the first patch covers the region close to the horizon while the second one is 
used in the asymptotic region (the same approach was employed in the recent work p4|]). 
The equations are relaxed on both patches one after the other. Moreover, the patches 
communicate, providing boundary data to each other in the overlapping region. 

We have considered the problem of constructing KK black holes by following the ap- 
proach in ||3^, [35|. A generalization to any D of the metric ansatz used there is given 
by 

ds^ = e2-§(-.-) J(r, z)Bl{r){dr'^ + dz^) + e^^^'^''^ R^{r, z)dujl_^ - e^^Al{r)dt^ , (A.7) 
where 

Ao(r-) = tanh {{D - 3)(r - ro)/2) , Bo(r) = (1 + e"(^-3){r-ro))2/(D-3)^ ^^ ^^ 

and 

J{r,z) = ^= f {AS 
2\/ sinh r -|- cos^ z 

R{r, z) = r + — log(l -|- 2e~'^ \/ sinh^ r -|- cos^ z -|- e''^'^ \j —\ -|- (e^'' -|- 2e^\/ sinh^ r -|- cos^ z)^ 
while ro < is an input parameter corresponding to the event horizon radius. As it was 



extensively discussed in [33|, tq also specifies the size of the black hole. Without any loss 
of generality, the compactification radius is fixed here to be L = vr. 

One can see from ( [A.7D that the metric functions A,B,C in the generic ansatz ( [4.1] ) 
contain some background terms {e.g. A{r,z) = A{r,z) -|-logj4o(r) etc.). The functions 
Ao{r), BQ{r) assure that, for small black holes, the solution is approximated near the 
horizon by the Schwarzschild-Tangerlini metric in D— dimensions. The background terms 
J{r,z), R{r,z) are used to implement the topology change between event horizon and 
spatial infinity. A detailed discussion of the choice ( |A.9| ) is presented in [^] . Here we only 
note that i?(r, 7r/2) = 0, for r < (thus also for r = tq), and R{r, Tr/2) > elsewhere, 
which shows that z = 7r/2 with r < represents the symmetry axis, while the topology 
of the event horizon is S^^"^. Also, J{r,z) diverges at r = 0, 2 = tt/2, which shows that 
e~'^^ there. However, this is just a coordinate singularity, as found e.g. by computing 
the Kretschmann scalar of the solutions. 

The numerics is performed with the functions A,B,C which are finite everywhere. 
The equations satisfied by these functions are derived again from the Einstein equations 
G* = 0, Gl + = and Gq = 0. As in the NUBS case, the remaining Einstein equations 



of a black ring is S^~^ x {i.e. it looks like a black string in the horizon region), the background 
being approached in the far field |^^. In the vacuum case, the only exact solution here corresponds to the 
D = 5 black ring in |5^. However, some of the black objects with a nonspherical topology of the horizon 
can be studied in a nonperturbative approach by employing a suitable version of the Weyl formalism 
IpOl]. This leads to a much simpler numerical problem as compared to the caged BHs in KK theory. 
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= 0, — = yield two constraints, which are automatically satisfied for the set of 
boundary conditions chosen. 

For the metric ansatz ( [A.7| ), the caged BH solutions are constructed subject to the 
following set of boundary conditions 

^1 =B\ =C\ =0, 

lr=oo lr=oo lr=oo ' 

drA\ =0, drB\ =1-^1 , d,C\ =1-^1 , (A.IO) 

' \r=ro ' ' lr=ro 9/ '''=''"0 Ir—r-Q D lr=ro' ^ ' 

^ \z=0,it/2 ^ l2=0,7r/2 b=0,7r/2 

The regularity of the solutions requires also B = C for r < 0,z = 7r/2. 

The Hawking temperature and entropy of the BH solutions are given by^^ 

Th = S=^f;dz e^(^0'^)+(^-3)f^(''"'^)j(ro,z)i?^-3(ro,z).(A.ll) 

The mass and tension of the solutions can be read in principle from the asymptotics of 
the metric functions A, B (in practice, these quantities are computed usually from the 



event horizon quantities, by using the Smarr relation (4.3) together with the first law of 
thermodynamics ( [4.4D ). 

The numerical iteration starts with large values of |ro|, i.e. small black holes, which 
are well approximated by the Schwarzschild-Tangerlini solution in D— dimensions. These 
configurations are used as initial data for the solutions with a smaller absolute value of tq. 
This results in a branch of caged black holes, which stops to exist for a critical value of 
ro < 0. 

We have constructed in this way a number of D = 6 solutions, with relatively small 
values of tq. In our approach, the numerical calculations were based on the Newton- 
Raphson method and are performed with help of the program FIDISOL [^], which provides 
also an error estimate for each unknown function. The numerical accuracy of the BHs was 
always lower than in the NUBS case, the errors increasing around the point r = 0,z = 7r/2. 
However, the results for the solutions we have constructed so far are in good agreement 
with those in [^] , [^] . Typical profiles of the metric functions A,B,C are shown in Figure 
4 (left). 

B. Slowly rotating, charged uniform black strings in Einstein-Maxwell 
theory 

For the parametrization considered in this paper, the uniform, charged black string solution 
in the EM theory is given by 

dsl = (r) (-^ + dz'^ + r'^ide + sin2 Odip'^ + cos^ 9diol_^)^ - ^dt"^, A = V{r)dt,(B.l) 
where 



^Note that the Einstein equation 0^—0 assures that e'^''^'^' ^^^'"^ / ^/jJr^Tj is constant on the horizon. 
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ro > being the event horizon radius and < U < 1. 

We consider a slowly rotating ansatz which follows from (|6.2|): 

ds'^ = dsQ — 2asin^ 9j{r)dtd(p, A = V{r)dt + ah{r) sin^ 6dip, 

where a is an infinitesimal parameter. The functions h{r),j{r) are solutions of a system of 



second order ordinary differential equations which are found from (6^) by replacing there 
^ = - logS7, B = logQ/{D -3), C = \ogn/{D -3). 
These equations admit the following exact solution 



Qe 



Jl^i ^'^n rD-^{u^{D - 3){D - 2)' ' U - 4' - 3 ' D - A ' - l\ ^ 
(I)_4)(D_1)^2_(^_5)P_2) /D-2 2D -A 2{D - 3) rp g_4 

y'2{D-2){D-3)U ' ' L»-3 ' D - A ' [/^ - 1 ^ r ^ 

where 2-?^i is the hypergeometric function. An unexpected feature here is the occurrence 
for D = 5 of a logarithmic term in the expression of J(r), 

= (l - '-if^^i^) - '-^^^$P^r^^'{r) + ^^r-303(r)l0gl7(r)) . (B.4) 

The angular momentum and magnetic dipole moment are read off from the asymptotic 
form of j(r) and h{r) (see (|6.7D , (|6.8D ) 



Qe(6 + 4[/2 + D(Z)-5)(l + ^^)) 1 Ur) = -_9£ 

^2{D -2){D -3){D -2)U r^-^ r^" 



The expression of the gyromagnetic ratio ( |6.10 ) derived in this way is (here we use also 



( 5.29 ) for the mass to charge ratio): 



2C/2 

^ = [D- 2){D -3) + {D-l){D- A)U^ ' ^^'^^ 



which varies in the range 



D-3<,<'°-^'''°-'". (B.6) 

^ D{D-b)+h ^ ' 



Thus, different from the corresponding case of (-asymptotically Minkowski) black holes in 
D-dimensions, the gyromagnetic ratio of the slowly rotating black strings in EM theory is 
not constant even in the first order of the perturbation theory. This can be understood 
as follows. After performing a KK reduction with respect to the Killing vector d/dz, the 
static EM uniform black strings (|B.1| ) become black hole solutions in a Einstein-Maxwell- 



dilaton theory in D — 1 dimensions, with a value a = 1/ ^J2{D — 2){D — 3) of the dilaton 
coupling constant. However, the results in |63] show that the inclusion of a dilaton field 
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in the action modifies the value of g for the slowly rotating black holes in (i— dimensions, 
which deviates from the standard result d — 2. This explains the relation ( |B.5D we have 
found for charged black strings in D— dimensions. 

We close this part by remarking that a similar study can be performed for slowly 
rotating charged solutions with all angular momenta degrees of freedom excited. The 
results we have found for a single angular momentum here (and in Section 6) are expected 
to hold there also. In the first order of perturbation theory, all gyromagnetic ratio would 
be equal and given by the relation ( p.5| ) in the UBS limit. 
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